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The asymptotic behaviour of the survival or reunion probability of vicious walks with short-range 
interactions is generally well studied. In many realistic processes, however, walks interact with a 
long ranged potential that decays in d dimensions with distance r as r"''"''. We employ methods of 
renormalized field theory to study the effect of such long range interactions. We calculate, for the 
first time, the exponents describing the decay of the survival probability for all values of parameters 
a and d to first order in the double expansion in e = 2 — d and 5 = 2 — d — cr. We show that there 
are several regions in the a — d plane corresponding to different scalings for survival and reunion 
probabilities. Furthermore, we calculate the leading logarithmic corrections for the first time. 
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INTRODUCTION 



Systems consisting of diffusing particles or random 
walks interacting by means of a long-range potential are 
non-equilibrium systems, which describe different phe- 
nomena in physics, chemistry and biology. From a phys- 
ical perspective they are used to study metastable super- 
cooled liquids [l|, Hj melting in type-II high-temperature 
superconductors [3|, electron transport in quasi-one- 
dimensional conductors [3| and carbon nanotubes @. 
From a chemical viewpoint the interest in these systems 
lies in the fact that some diffusion-controlled reactions 
processes rely on the diffusion of long-range interacting 
particles which react after they are closer than an effec- 
tive capture distance. Some examples include radiolysis 
m liquids [6^, electronic energy transfer reactions [7| and 
a large variety of chemical reactions in amorphous me- 
dia [8|. From a biological viewpoint, the investigation of 
these systems is helpful in understanding the dynamics of 
interacting populations in terms of predator-prey models 
[^■[To| and membrane inclusions with curvature-mediated 
interactions pll [l^. 

Vicious walks (VW) are a class of non-intersecting ran- 
dom walks, where the process is terminated upon the 
first encounter between walkers p^ . The fundamental 
physical quantity describing VW is the survival proba- 
bility which is defined as the probability that no pair 
of particles has collided up to time t. Diffusing parti- 
cles or walks that are not allowed to meet each other 
but otherwise remain free, we call pure VW. The be- 
havior of pure VW is generally well-known. The sur- 
vival probability for such a system has been computed 
in the framework of renormalization group theory in ar- 
bitrary spatial dimensions up to two-loop order [I4l4l6l |. 
These approximations have been confirmed by exact re- 
sults available in one dimension from the solution of the 
boundary problem of the Fokker-Plank equation [l^] , 
using matrix model formalism [l3| and Bethe ansatz tech- 
nique [l^. On the other hand the effect of long range 
interactions has been extensively investigated in many- 
body problems. It has been shown that the existence of 
long-range disorder leads to a rich phase diagram with 



interesting crossover effects [I9l - l21| . If the potential is 
Coulomb-like (~ r~^~"') then systems in one dimension 
behave similar to a one-dimensional version of a Wigner 
crystal |2d| for ct < and similar to a Luttinger liquid 
for cr > [23|. If the potential is logarithmic then in the 
long-time limit the dynamics of particles are described 
by non- intersecting paths [13, HI] ■ The generalization of 
VW that includes the effect of long range interactions 
has not attracted much attention in the literature. Up 
to our knowledge there was one attempt to study long- 
range VW [IHl. Here the authors considered the case of 
a long-range potential decaying as gr~°'~'^, where g is a 
coupling constant. It was shown by applying the Wilson 
momentum shell renormalization group that only one of 
the critical exponents characterize long-range VW. For a 
specific value of cr {a = 2 — d) they show that the ex- 
ponent 7, which determines the decay of the asymptotic 
survival probability with time, is given by the expression: 



p{p~ 1) 
7 = z "1, 



(1) 



where p the number of VW in the system, ui — (e/2 + 
[{e/2)'^+gY/^) and e = 2 — d. There are limitations to the 
above approach. First, it is restricted to a single form of 



the potential 



') and systems such as membrane in- 



clusions and chemical reactions have different power-law 
potentials. Second, it considers identical walkers but one 
would like to have results if the diffusion constant of all 
walkers are different. Finally it is not convenient to com- 
pute higher-loop corrections using the Wilson formalism. 

In this paper we reconsider the problem of long-range 
VW using methods of Callan-Symanzyk renormalized 
field theory in conjunction with an expansion in e = 2 — d 
and S = 2 — d — a. We note that it is more convenient to 
compute logarithmic and higher loop corrections by using 
this method. We derive the asymptotics of the survival 
and reunion probability for all values of the parameters 
(ct, d) for the first time. 

In this paper we will show that there are several re- 
gions in cr — d plane in which we have different behavior 
of the critical exponent. Our results are summarized in 
Table I. We note that results on the line a + d — 2 have 
been obtained before [25|. Regions I and IV correspond 
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TABLE I: One-loop survival probability of p sets of particles 
with Tij particles in each set large-time asymptotic at different 
regions of the a — d plane. We refer to Figure 3 for specific 
value of a and d in each region. 



Region 



Survival probability 



I 

II 
III 
IV 

V, d = 2 
VI, a = 2 - 



^-(d-2)/2 _|_ ^-(d+<T-2)/2 

^-(d-2)/2 
i-^E„"i"j(l+«/21ogt) 



"Ml is defined by the formula ITl . 



to Gaussian or mean- field behavior (see Figure 3). In re- 
gion II we found that the system reproduces pure VW. 
Logarithmic corrections in region III and at the short- 
range upper critical dimension d — 2 have been obtained 
as series expansion in 5 = 2 — cr — d. 

The remainder of this paper is organized as follows: 
Section |lT] reviews the field theoretic formulation of long 
range VW and describes Feynman rules and dimension- 
alities of various quantities. In section Hill we derive the 
value of all fixed points and study their stability. Sec- 
tion IIVI presents results for the critical exponents and 
logarithmic corrections of various dynamical observables. 
Section|V]contains our concluding remarks. In Appendix 
A we give the details of the computation of some integrals 
that appear in Section IIIII 



II. MODELLING VW WITH LONG-RANGE 
INTERATIONS 

As the starting point of the description of our model we 
consider p sets of diffusing particles or random walks with 
Ui particles in each set i = 1 . . . p, with a pairwise intraset 
interaction which includes a local or short-range part and 
a non-local or long-range tail. The local part determines 
the vicious nature of walks: if two walks belonging to 
the different sets are brought close to each other, both 
are annihilated. Walks belonging to the same set are 
supposed to be independent. At t = all particles start 
in the vicinity of the origin. We are interested in the 
survival and reunion probabilities of walks at time i > 0. 

A continuum description of a system of N Brownian 
particles Xi with two-body interactions is simplified by 
the coarse-graining procedure in which a large number of 
microscopic degrees of freedom are averaged out. Their 
influence is simply modelled as a Gaussian noise-term 
in the Langevin equations. A convenient starting point 
for the description of the stochastic dynamics is the path- 



integral formalism. Then the system under consideration 
is modeled by the classical action 



S= dtK]Af/(2A) + E^(^^-^^- 



(2) 



where t is (imaginary-)time, Xi{t) is the d-dimensional 
vector denoting the position of ith particle at time t. 
Di is an ith particle diffusion coefficient. The path- 
integral representation of the probability density function 
for the particle displacements from their original posi- 
tions is given by the functional Z — J T>X exp[—S]. The 
survival probability is defined as the expectation value 



pit)^{Y[[i^6{x,{t)-x,m) 



(3) 



with respect to the functional Z. It is computed in the 
framework of usual perturbation theory and will be a sum 
of integrals over internal degrees of freedom. It is more 
convenient to perform these integrations in Fourier space. 
To do this we would need the Fourier transform of the 
interaction potential V{r). We note that it is comprised 
of a short-range part of the form Vb(r) = XS{r) and a 
long-range part which decays with the distance r as a 
power law, Vi{r) = gr^'^^'^ . The Fourier transform of 
the latter is divergent if ct > 0. We introduce the cut-off 
parameter a to regularize the singularity Vi{r) — g(r^ + 
^2^-(d+(j)/2 ^ Fourier transformation of this function is 
given by the expression 

= g^^^iQ/^r^'K^Maq), (4) 

where Kg. is the modified Bessel-function with index a. 
Small a expansion of (|3]) at leading order yields 



log(aq). 



if 
if 



cr = 



(5) 



where we used the property K-a-ix) — K„{x) of the 
Bessel function. The non-universal coefficient coming 
from the Taylor expansion can be absorbed by the ap- 
propriate renormalization of the constant g. Special cases 
when a is even gives logarithmic behavior. Effectively it 
does not change our results. So we focus on the typical 
term . 

The second quantized version of the action ^ can be 
constructed using standard methods [1^ [l^l . The gener- 
alization of the action to the long-range interacting case 
is also known [29, . The result is 
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5((/)„</)t) = / dtd''x{Y,[4^tc|>^+D,V<|>lWc|,,]}+ / dtd 




FIG. 1: Feynman rules for the theory (|6|. Notice that both A 
and g vertices appear with differenet i and j indices and that 
g has momentum dependence. 

The first term describes the evolution of free randora 
walks vi^ith diffusion constants Di. The potential is 

Vijix - y) = XijS{x ~y)+ gtjV{x - y), (7) 

and we refer to Xij , gij as short-range and long-range cou- 
pling constants respectively. 

A dynamic response functional associated with the ac- 
tion ([6]) is 

Z = J r>,^I?(/ite-^('^-'^.) (8) 

where 0^ {x, t) is the complex scalar field. After the 
quantization we may treat (x, t) as the creation op- 
erator which creates a particle of sort i at point x at 
time t. Having the dynamic response functional, corre- 
lation functions can be computed as functional averages 
(path integrals) of monomials of and 0^ with the weight 
exp{-5(0,0t)}. 

As a first step towards the renormalization group anal- 
ysis of this model, we discuss the dimensions of various 
quantities in ^ expressed in terms of momentum: 

[t]=p-^ [<^]=/ [A]-p2-<i [g]^p^-d-. (9) 

The naive dimension of the coupling constant g allows us 
to identify the upper critical dimension dc{(T) = 2 — d—a. 
For cr > 0, the short-range term naively dominates the 
long-range term and we expect to have the behavior of 
the system similar to the case of pure VW. We will re- 
serve the symbol e (e = 2 — d) to denote deviations 
from the short-range critical dimension dc = 2, and 5 



xd''yY,cl^lit,x)Mt,x)V,jix - y)^]{t,y)cj),{t,y) . (6) 

i<j 

I 

(S = 2 — d—ij) for the deviations from the long-range crit- 
ical dimension dc{cr). If u = then the critical dimension 
of the long-range part coincides with the short-range part 
and we have the non-trivial correction to the asymptotic 
behavior due to long-range interactions. This bound- 
ary separates mean-field or Gaussian behavior from long- 
range behavior. For a < the long-range term dominates 
the short-range term and we expect to have non-trivial 
corrections to the behavior of the system. 

Now we consider diagrammatic representation ele- 
ments of model In zero- loop approximation the ver- 
tex 4-point function takes a simpler form after Laplace- 
Fourier transformation: 

r[f'\s,p) = v,,{p^+p2)S{Y,p,). (10) 

k 

The same transformation applied to the bare propagator 
yields: 

Tf^'\s,p) = {s + D.,p'r^ (11) 

We note that there are no vertices in ^ that produce di- 
agrams which dress the propagator, implying there is no 
field renormalization. As a consequence the bare propa- 
gator pT|) is the full propagator for the theory. Feynman 
rules are summarized in Figure I. There are two ver- 
tices in the theory: one is a short-range A-vertex and 
another is a long-range momentum dependent g-vertex. 
Each external line of the vertex corresponds to a func- 
tionally independent field. The propagator is formed by 
contracting appropriate lines from different vertices. We 
recall the propagator is the correlation function of 0i and 
fields only. 

Physical observables are computed with the help of 
correlation functions. The probability that p sets of par- 
ticles with Hi particles in each set start at the proximity 
of the origin and finish at Xi^ai (* index enumerates dif- 
ferent sets and index enumerates particles in set i) 
without intersecting each other can be obtained by gen- 
eralizing eqn ([3]). In the field theoretical formulation, this 
probability becomes the following correlation function: 

n n rf'a;,,„,(0,(i,a;.,„J(0l(O,O))"'), (12) 

In the Feynman representation it is the vertex with 2A^ 
{N ~ J2j ^i) external lines. In the first order of the per- 
turbation theory one needs to contract these lines with 
corresponding lines of the vertices in Figure 1. Since 
there are many independent fields in the correlation func- 
tion (jl2[) this operation can be done in many ways. It 
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FIG. 2: One-loop Feynman diagrams contributing to Xmj 



yields a combinatorial factor, ninj, in front of each di- 
agram, which is the number of ways of constructing a 
loop from the rii lines of type i and rij lines of type j on 
the one hand and one line of type i and one line of type 
j on the other hand. From the next section we will see 
that the survival probability scales as G{t) ~ t~'' , where 
7 is the critical exponent. If all walks are free, 7 = 0. 
In the presence of interactions we expect 7 to be a uni- 
versal quantity that does not depend on the intensity of 
the short-range interaction Xij . It is convenient to intro- 
duce the so called truncated correlation function which 
is obtained from by factoring out external lines: 



(13) 



Another physical observable, the reunion probability, 
is defined as the probability that p sets of particles with 
rii particles in each set start at the proximity of the ori- 
gin and without colliding into each other finish at the 
proximity of some point at time t: 



convenient scheme developed for dealing with these diver- 
gences follows Callan-Symanzik renormalization-group 
analysis ^,'31^. Within this scheme we start with the 
bare correlation function G{t] A, g), where A = {A^ }, and 
9 — {9ij} denote the set of bare short-range and long- 
range coupling constants. In the renormalized theory it 
becomes Gi^(t; A^j, g^j, /x). From dimensional analysis it 
follows that 



(15) 



where fj, is the renormalization scale. The scale invariance 
leads to the expression 



Gnit; Xr, gii, pi) = Z{Xii, g^, ii)G{t; A, g). 



(16) 



Here functions Z are chosen in such a way that 
Gnit, Xr, gR, I) remains finite when the cut-off is removed 
at each order in a series expansion of Xr, g^, e and 
6. From the fact that G{t, A, g) does not depend on 
the renormaUzation scale /x we get the Callan-Symanzik 
equation 



d_ 

9/i 



In the Feynman representation it is depicted as the wa- 
termelon diagram with 2N stripes. We note that if the 
theory is free this expression is the product of free prop- 
agators and at the large-time limit the return probability 
scales as Ro{t) ^ i-{N-i)d/2 ^ interactions are taken 
into account it becomes R{t) ^ f~{N-i)d/2-2j ^ where 7 is 
survival probability exponent. The reason that it enters 
with the factor 2 is the following. If we cut a watermelon 
diagram of the reunion probability correlation function 
in the middle then it produces two vertex diagrams with 
2N external lines of the survival probability correlation 
function. As a result the reunion probability is the prod- 
uct of two survival probabilities. It remains true in all 
orders ofperturbation theory. For a rigorous proof we 
refer to [191 . 




III. THE RENORMALIZATION OF 
OBSERVABLES 

While computing correlation functions like pertur- 
batively one faces divergent integrals when d — dc- The 



FIG. 3: The critical behavior of vicious walks with long-range 
interactions in the different regions of the {a, d) plane. Region 
I and IV correspond to the mean field short-range behavior, 
in region II will be critical short-range behavior, region III 
is the long-range behavior. The Unes d — 2 and a + d — 2 
represent regions V and VI respectively. 
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where the /3-functions are defined by 



d 

and the function 7 by 



(18) 



(19) 



The renormahzation group functions are understood as 
the expansion in double series of couphng constants A 
and g and deviations from the critical dimension e and S. 
We take S = 0{e). The coefficient Z{Xji, g^, ^) is fixed 
by the normalization conditions. It is more convenient 
to impose these conditions on the Laplace transform of 
the truncated correlation function (|13l) . One sets the 
following condition then 



(20) 



when s = fj,. We note that the same multiplicative renor- 
mahzation factor Z yields T finite. From this fact one 
can infer that 



r(Ai;A,5) = z{m\,gy 



(21) 



If we express unrenormalized couplings in terms of renor- 
malized ones Ij2ip we will obtain the equation for finding 
Z explicitly. 

The equation (fT7|) can be solved by the method of char- 
acteristics. Within this method we let couplings depend 
on the scale which is parametrized by /i(a;) = x/x. Here 
X is introduced as a parametrization variable of the RG 
flow and is not to be confused with position. Henceforth 
X will refer to this parametrization variable. We intro- 
duce running couplings X{x) and g{x). They satisfy the 
equations 

a;^5(a;) = Pg{X{x),g{x)) x-^X{x) = l3x{X{x), g{x)). 

. (22) 

The renormalized value should be defined by the initial 
conditions A(l) = Xji and g{l) = g^. the solution of the 
equation is then 

GR{t)^e^ GR{^l-^;X{^lt),g{^lt),fl) (23) 

Next we calculate the first-order contribution to the 
renormalized vertices. The A-vertex is renormalized by 
the set of diagrams that are shown in Figure 2. We notice 
that there are no diagrams producing the momentum de- 
pendent (7-vertex in the theory (jG)). This statement is the 
corollary of the fact that only independent fields of power 
one enter into the expression of the vertex and there are 
no higher powers of fields. Also we keep in mind that 
the renormalized couplings are defined by the value of 
the vertex function taken at zero external momenta. It 
produces the following expression: 



XRij 

gRij 



= A 



^3 

9tj 



\{Xlh+2\jg,jh+glh 



(24) 



where Ik = Ik{o", Di, Dj) are one-loop integrals corre- 
sponding to the diagrams a, b, c in the Figure 2 respec- 
tively. Using the Feynman rules we can explicitly write 
them down: 



(fc-l)cr 



fc = l,2,3. (25) 



We will use dimensional regularization procedure to com- 
pute these integrals. The details of the computation are 
summarized in Appendix A. We note that integrals will 
diverge logarithmically at different values of the spatial 
dimension d. For this reason it leads to different critical 
behavior in different regions of the a — d plane (see Fig- 
ure 3). These regions correspond to four possibilities for 
e = 2 — d and ^ = 2 — d — cr to be positive or negative. 
Only if (5 = 0(e) or, in other words, if both e and S are 
infinitesimally small but the ratio e/S is finite we expect 
non-zero fixed points of the renormahzation group fiow. 
Similar approximation have been used before 19] but for 
different models with long-range disorder. It allows us to 
follow the standard procedure of deriving the /3-functions 
which consists of two steps. 

First, we express unrenormalized couplings in terms of 
the renormalized. For the short-range coupling constant 
A it can be done by solving the quadratic equation in 
(1241) . Expanding the square root and keeping terms up 
to the second order we infer that 



I Xij 
\9ij 



^Rij 

gRi] 



2 y^Rij ' 



"^^RijgRij^ 



9Ri] 2.5-e) 

(26) 

where a^, bd and Crfcoefficients have been found explicitly 
in Appendix A. Now we introduce dimensionless renor- 
malized couplings 



gRtj = ad{2s) ■'/^ XRij = bd[2s) 



-ell 



(27) 



An important observation is that Cdad = b^ which can be 
verified by explicit substitution (see Appendix A). Mul- 
tiplying the first and second equation in (1261) by the fac- 
tors ad and bd respectively, and using redefinitions (j27p 
we can condense all pre-factors in the right hand side of 
the equations into the dimensionless constants. 

Second, we differentiate equations (I^Bl) with respect to 
the scaling parameter /i. Using definitions (|18p and the 
fact that bare couplings do not depend on the scale, we 
derive 




= —£^R.ij + (Aflij + gRijY 

= -SgRij 



(28) 



where the right hand side is understood as the lead- 
ing contribution to the /3-functions from the double ex- 
pansions in A, 5 and £,(5. From (j28p we see that it is 
convenient to introduce new coupling constants URij ~ 
^Rij + gRij- After this step the renormahzation group 
equations read 
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f3g,ij = -SgRij 

We note that in the last equations g couphng constant 
has been redefined agmj — ?► gmj- 

Fixed points are zeros of the /3-functions. If (5 ^ then 
the last equation in (|29p is zero only when g^, = 0. Then 
the first equation has two solutions u = and u ^ e. 
If (5 = then g plays the role of a parameter and the 
fixed points are determined by the roots of the quadratic 
equation 

= -eu + - g (30) 
which are real if <? > — (e/2)^ and we find 



IV. CALCULATION OF CRITICAL 
EXPONENTS AND DISCUSSION 

Here we describe our method of computing critical ex- 
ponents. It is based on the formula from the previ- 
ous section. First, we obtain the leading divergent part 
of the correlation function. The renormalized correlation 
function depends on the scale /i but it appears in all for- 
mulas in combination with time: fit. Second, since we 
have found the bare coupling constant as a function of 
renormalized (dressed) couplings we express correlation 
function in terms of dressed couplings. Finally using the 
normalization condition ([^0]) and the definition ([T^ we 
differentiate Z with respect to fid/dfj, to obtain the ex- 
ponent 7. The poles should cancel after this operation. 

In section 2 it was explained that the truncated cor- 
relation function in the one-loop approximation is given 
by the formula 



ui,2 = e/2± V(e/2)2+.g. 



(31) 



All fixed points are listed in the Table II. The stability of 
these fixed points is determined by the matrix of partial 
derivatives 



^ __{ dpu/du dPuldg \ 
P* [ dpjdu dPg/dg j 



(32) 



u=u,,g=g. 



Eigenvalues are listed in the Table HIl The Gaussian fixed 
point is stable in all directions for e < and 5 < which 
corresponds to region I in Figure 3. In this region we 
find both short-range(pure VW) and long-range mean- 
field behavior depending on the sign of a. On the con- 
trary, for e > and (5 > we find that the Gaussian 
fixed point is unstable(irrelevant) in all directions and 
the short-range (pure VW) fixed point is stable(relevant) 
only in u-direction. It means that long-range interactions 
will play a leading role. This region corresponds to re- 
gion III in Figure 3. Next for e > and (5 < we find 
that the short-range (pure VW) fixed point is stable in 
all directions. It means that the system is insensitive to 
the long-range tail. This region corresponds to region II 
in Figure 3. Finally for e < and (5 > we find that 
the short-range (pure VW) fixed point is unstable in all 
directions and the system will be described by mean-field 
at long time. 



TABLE II: Fixed points for flow equations (|29|l and the cor- 
responding eigenvalues (Ai,A2) of the stability matrix (|32|l . 
We note that u\ and U2 are values of the 



Fixed point 


(u.,g.) 


(Ai,A2) 


Gaussian 


(0,0) 


(e,5) 


Pure VW 


(e,0) 


{-e,5) 


LR stable 


(ui,0) 




LR unstable 


(^^2,0) 


(Ve2-4»,0) 



T{t\ A, 5) = 1 - ^ UiTij [Xijh + gijh 



(33) 



Here integrals are the same as in d^S]) . 

We start our analysis with the region I. Notice that 
truncated correlation function T{t) and survival prob- 
ability G{t) have similar large time behavior. We use 
large momentum cut-off to compute integrals Ii and I2 
as in formula (I64[) in Appendix A. The renormalization 
of coupling constants is trivial in this case. Therefore the 
leading contribution to the survival probability is given 

by 



G{t)^t^^^-^^/^+got^^-''-'''i/\ 



(34) 



where g^ is non-universal coefficient and we will not need 
its exact value. We notice that if cr > the second term 
will decay faster than the first term and in the long-time 
limit it will produce the same behavior as mean-field pure 
VW. On the other hand if ct < the first term will de- 
cay faster and long-range interactions will play a leading 
role. Many authors observed similar behavior in vari- 
ous systems with long-range defects [l9l - [2l| . Intuitively 
if potential falls fast with distance than the system effec- 
tively represent system with short-range potential where 
particle interact when they are close to each other. 

Region IV exhibits similar behavior. Now the integral 
I2 is computed with the help of the dimensional regular- 
ization (I55|) and the integral Ii remains the same. From 
the fact ([T5|) one can infer that the survival probability 
scales as 



(35) 



Short-range behavior dominates because the running 
coupling constant will flow towards the Gaussian fixed 
point at long time limit which is the only stable fixed in 
this region. This result is exact regardless the number of 
loops one takes into account. 
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In Region II the computation is as follows. 

In Z = ^ n,n, (a., ^ + g,,t^^-^-'^^'^) , (36) 

so plugging the result from ([59| to ([36| we obtain at the 
fixed point (A* =£,3 = 0) 



7 = 



(37) 



And we reproduce the pure VW behavior. This result is 
the reflection of the fact that the renormalization-group 
trajectories run away to stable pure VW fixed point. It 
is with agreement with the results obtained by Katori in 
[l7| for d — 1, and the logarithmic intraset particle inter- 
actions. The irrelevance of the long-range interaction in 
lower dimensions is a typical phenomenon observed in a 
various out of equilibrium interacting particle systems. 

We now consider regions III, V and VI. Integrals in ((33|) 
are computed via dimensional regularization. Taking the 
inverse of (j33l) and then logarithm one can obtain at the 
leading order: 



\ogZ = niUj f Ay — -I- gij 



(38) 



where and bd are defined in Appendix A in (j59p and 
([60| . We note that after taking the derivative the poles 
in (1551) will cancel in the limit of S — 0{e). Also one 
recalls the expansion (1^^ and the redefinitions in (P7| . 
Using pQ]) we show that the expression for the function 
7 which determines critical exponent takes the form 
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niUjUR 



(39) 



Evaluated at the stable fixed point (ui 
a/ (e/2)^ -I- g it gives the following result: 



e/2 



(40) 



and the survival probability scales as G{t) ~ 

We will now find the logarithmic corrections to this 
scaling law. The running coupling constant can be found 
from the flow equation (P^ : g{x) = e~^^g. In the case 
5,e = (the intersection of regions V and VI) the flow 
equation for u{x) is 



du(x) 
dx 



-u^{x) + g 



(41) 



and the solution is 



u{x) = y/gta.nh{y/g\ogx + (j)o) y/gtanh{y/g\ogx), 

(42) 

where 0o is the initial condition and we do not need its 
exact form. After plugging this expression into the (p3)) 
we infer 

fit 

[ - - dx 

j l{u,g) — ^ log(cosh(y^log^<)) (43) 



Thus the survival probability is 

G(0^cosh(vglogt)-^^' 



(44) 



In the limit of large time cosh(y^logt) ^ imply- 
ing gamma = —^J2ij^i''^jV9 which is consistent with 
equation (HUl) . For negative coupling constant g < the 
solution in (l42l) becomes 



(45) 



u{x) 



g|tan(V|g|loga;) 



The integral (j43| is divergent if t > exp(7r/2-\/[gf) which 
leads to the result that the survival probability is zero 
beyond this time. For smaller times one has G{t) ^ 
cos(-\/|5|logt)~ 2 5Z . Thus, upto one- loop order ap- 
proximation. It implies that if walks are attracted to each 
other then all of them will annihilate at some finite time. 
This might be a signature of faster than power law de- 
cay and we expect to have corrections to this behavior at 
higher loop approximation. 

Next we consider the case when e — and S ^ but 
6 remains small i.e. region V. The flow equation for the 
u{x) is 



X — ; — = —u (X) + gx 
dx 



(46) 



and the solution can be found by the method of pertur- 
bation. Up to the flrst order 

u{x) = ^tanh(^loga;) + (5V5log(a^) tanh(^logx) 

(47) 

After plugging this expression into eqn ()23p we infer 



iiu,g) 



dx 



(48) 

Therefore we have the correction to the survival proba- 
bility in the form 



Q ^ ^-|E«i»i\/9(l+<5/2(logt)) 



(49) 



Now we extend our analysis to the case when e > 0, 
corresponding to regions III and VI. The evolution of the 
coupling constant is 



x—u(x) = eu — + qx^ 
dx ' 



(50) 



We choose the ansatz in the form u{x) — uq{x) + dv{x). 
For (5 = (i.e. region VI) the equation for ^0(2;) reads 



d , ^ 2 
x—uq(x) = ewo - + 5 



(51) 



and we reproduce the result P0|) . We now extend to the 
case where e,5 > (region III). Here we will need the 
exact solution to ([5T|) to flnd the corrections: 



Uo{x) 



(52) 
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where C = {un — U2)/{ui — un). The logarithmic cor- 
rection follows from the form of the perturbation. The 
equation for v{x) is 

x—v{x) = ev ~ 2uqv ~ glogx (53) 

The solution can be found explicitly as a combination of 
hypergeometric functions. In the most interesting case, 
e = 1 {d = 1) the hypergeometric functions are degen- 
erate and become linear functions. Corrections to the 
integral then read 

Jjdx/x ^ ^Su,\og\t)+\og{t){tr~"-) (54) 
1 

In the limit of large time only the first term contributes 
to the exponent and the survival probability scales as 



evaluated in [13] by using alpha representation l/((j^ + 
s) = J^°° dae^^'^ and the result is 

h - Kd{2s)-''^T{e/2). (56) 

We notice that since there is no angular dependence one 
can perform d ~ 1 integrations and one will be left with 
one dimensional integral. To compute this integral we 
use the formula [3^ : 



We see that in our case P = s, Q — [Di + Dj) and 
h' ~ d + (ii — 1)(T. This immediately gives the result: 



CONCLUSION 



(55) 



In summary, we studied long-range vicious walks using 
the methods of Callan-Symanzik renormalized field the- 
ory. Our work confirms the previously known RG fixed 
point structure including their stability regions. We cal- 
culated the critical exponents for all values of a and d to 
first order in e expansion and to all orders in i5 expan- 
sion, which have hitherto been known only for d + a = 2. 
Our results indicate that, depending on the exact values 
of d and cr, the system can be dominated by either short 
range (pure VW) or long range behaviors. In addition, 
we calculated the leading logarithmic corrections for sev- 
eral dynamical observables that are typically measured 
in simulations. 

We hope that our work stimulates further interest in 
long-range vicious walks. It would be interesting to see 
further simulation results for the critical exponents for 
d > 1 and for logarithmic corrections. Also, it would be 
interesting to have analytical and numerical results for 
other universal quantities such as scaling functions and 
amplitudes. 



VI. ACKNOWLEDGMENTS 

AG would like to acknowledge UC Merced start-up 
funds and a James S. McDonnell Foundation Award for 
Studying Complex Systems. 



Appendix A 

Effective four-point function (one-particle irreducible, 
IPI) that appeared in (IM)) is composed of usual short- 
range and new momentum dependent vertices. This gives 
rise to integrals The first integral /i = 1 has been 



1 



S 2 



2 \{D,+Dj 



where Ka = 2'^-^TT-'^^^r-^{d/2) is the surface area of 
d-dimensional unit sphere. 
It is convenient to define 



ad = 



2 \{D,+D,) 



d/2 



{2s) 



-e/2 



bd^-^ 



2 V(A + ^j) 



{d+a)/2 



i2s) 



-S/2 



Kd 



Cd 



(d+2(j)/2 



(2s) 



-{25-e)/2 



(59) 



(60) 



(61) 



2 \{D^ + Dj), 
So integral in the limit oi d — 0{s) can be written as: 

We used an expansion r(e/2) ^ 2/e for small e. An 
important property of coefficients ([5^)) - (pT|) is that 



which can be verified by direct substitution. 
Now we compute mean field integrals: 



(63) 



= y d''qdtq''+'' exp(-i(A + D,)q^) ^ t-^d+a-2)/2^ 

(64) 

where we assumed that the large momentum cut-off is 
imposed and corresponding coupling constants have been 
renormalized. The non-universal coefficient is not impor- 
tant. 
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